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ABSTRACT 

Luminous X-ray outbursts with variability amplitudes as high as ~ 1 000 have 
been detected from a small number of galactic nuclei. These events are likely 
associated with transient fueling of nuclear supermassive black holes. In this 
paper, we constrain X-ray outbursts with harder spectra, higher redshifts, and 
lower luminosities than have been studied previously. We performed a systematic 
survey of 24 668 optical galaxies in the Chandra Deep Fields to search for such 
X-ray outbursts; the median redshift of these galaxies is ~ 0.8. The survey spans 
798 days for the Chandra Deep Field-North, and 1 828 days for the Chandra Deep 
Field-South. No outbursts were found, and thus we set upper limits on the rate 
of such events in the Universe, which depend upon the adopted outburst X-ray 
luminosity. For an outburst with X-ray luminosity > 10 43 ergs s _1 and a duration 
of 6 months, the upper limit on its event rate is ~ 10~ 4 galaxy -1 yr _1 , roughly 
consistent with theoretical predictions. Compared to previous survey results, our 
harder-band and deeper survey suggests that the outburst rate may increase by 
a maximum factor of 10 when considering both obscured X-ray outbursts and 
redshift evolution from z ~ to z ~ 0.8. Our results also suggest that the 
X-ray luminosity function for moderate-luminosity active galactic nuclei is not 
primarily due to stellar tidal disruptions. 

Subject headings: galaxies: active — galaxies: nuclei — X-rays 



1. INTRODUCTION 

X-ray observations, mainly with ROSAT, have found about six transient, large-amplitude 
X-ray outbursts from galactic nuclei (e.g., Donley et al. 2002; Komossa 2002; Komossa et al. 
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2004; Vaughan, Edelson, & Warwick 2004; and references therein). These events have vari- 
ability factors of ~ 20-1 000, peak 0.5 — 2 keV X-ray luminosities comparable to those of local 
Seyfert galaxies (~ 10 42 -10 44 erg s _1 ), decay timescales of months, and soft X-ray spectra. 
They have been observed in both active and inactive galaxies (i.e., galaxies with and without 
a persistently accreting supermassive black hole, respectively). The estimated event rates of 
these outbursts suffer from large systematic and statistical uncertainties. In nearby inactive 
galaxies the event rate is ~ 10~ 5 galaxy -1 yr _1 , while in active galaxies the event rate ap- 
pears to be ~ 100 times higher (e.g., Donley et al. 2002). At least some of these outbursts 
induce accompanying optical nuclear variability (e.g., Brandt, Pounds, & Fink 1995; Grupe 
et al. 1995). Recently, five new candidate X-ray outbursts were found in the XMM-Newton 
Slew Survey (jEsquej et al.l 120071 ). and one outburst has been detected in t he ultraviolet and 
is thought to have the same origin as X-ray outbursts (IGezari et al.ll2006l ). 



The physical origin of nuclear X-ray outbursts remains mysterious. The most likely 
explanation is that they are caused by transient fueling events of nuclear supermassive black 
holes (SMBHs). Since SMBHs are thought to be ubiquitous in nucleated galaxies (e.g., Fer- 
rarese & Ford 2005 and references therein), fueling events seem inevitable in crowded galactic 
centers. Fueling may occur when a star, planet, or gas cloud is tidally disrupted and partially 
accreted (e.g., Rees 1990; Ulmer 1999 ; Li, Narayan, &: Menou 2002). The predicted ey ent 



rate for stellar tidal disruptions (e.g., iMagorrian fc Tremaineill999l ; ISyer fc Ulmerlll999l ) is 



roughly consistent with the poorly constrained rate of large-amplitude X-ray outbursts in 
inactive galaxies. However, I Wang Sz Merritt J2004T ) predicted a rate that was a factor of 
~ 10 higher than earlier results, mainly due to downwardly revised black-hole masses. This 
predicted rate, apparently in excess of the observed rate of X-ray outbursts, may indicate 
that a large fraction of tidal- disrupt ion events do not exhibit expected X-ray outburst charac- 
teristics, perhaps due to short durations of the events or X-ray obscuration. In some cases, 
transient fueling could be due to accretion- disk instabilities (e.g., Siemiginowska, Czerny, 
& Kostyunin 1996). Some outbursts could also perhaps be explained by X-ray afterglows 
of gamma-ray bursts (GRBs), though no sim ultaneous GRBs were reported for the known 
X-ray outbursts (e.g., 



Komossa fc Baddll999h 



Aside from the innate scientific interest of nuclear X-ray outbursts, determining the rate 
and properties of such events is relevant to planning future missions such as the Black Hole 
Finder Probe (BHFP), Lobster, and eROSITA. The BHFP, for example, is likely to observe 
in the hard X-ray band, and the large sky coverage of this mission will allow intrinsically rare 
transient events to be studied (e.g., Grindlay 2005). Facilities such as the Large Synoptic 
Survey Telescope, the Joint Dark Energy Mission, and the Laser Interferometer Space An- 
tenna should also allow the accompanying optical (e.g., Brandt 2005) and gravitational- wave 
(e.g., Kobayashi et al. 2004) outbursts to be studied. 
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Previous studies of X-ray outbursts from galactic nuclei have delivered fascinating results 
but require extension so that this phenomenon can be better understood. For example, 
advances are required in the following directions: 

1. Outbursts with harder X-ray spectra. As noted above, the outbursts discovered to 
date have generally had soft X-ray spectra with effective 0.2-2.0 keV power-law photon 
indices of T m 3-5. However, this may partially be a selection effect owing to the soft 
X-ray bandpass of the ROSAT satellite. Transient fueling events of SMBHs should be 
capable of generating harder X-ray spectra (r 1.7-2.2) via Compton upscattering in 
the accretion flow; such spectra are observed from most active galactic nuclei (AGNs) 
as well as from transient Galactic black holes. Indeed, one X-ray outburst has shown 
evidence for spectral hardening as it declined (Komossa et al. 2004). In addition, 
obscuration can harden the observed X-ray spectra of an outburst, which would make 
it more difficult to detect in the current outburst surveys. A fairly small column density 
of Ah = 5 x 10 21 cm~ 2 reduces the expected 0.2-2.0 keV flux from a soft-spectrum 
outburst by a factor of ~ 20. 

2. Higher redshift outbursts. The X-ray outbursts discovered to date are at low redshift 
with z = 0.01-0.15. Since only a small fraction of cosmic time is spanned by this 
redshift range and the source statistics are limited, little is known about the redshift 
evolution of their frequency. It is plausible that the X-ray outburst rate could show 
evolution over cosmic time, considering that AGNs and galaxies both show strong 
evolution. The number densities of comparably X-ray luminous AGNs evolve upward 
by factors of 10-30 out to z ~ 1 (e.g., Brandt & Hasinger 2005 and references 
therein), and Milosavljevic, Merritt, & Ho (2006) have suggested that a significant 
portion of the X-ray luminosity function of such AGNs may be comprised of sources 
powered by tidal disruptions. 

3. Lower luminosity outbursts. Transient phenomena in complex natural systems often 
follow power-law distributions in frequency of occurrence, such that low-amplitude 
events are more common than high-amplitude events (e.g., earthquakes, avalanches, 
solar flares, and ecological extinctions). It is plausible that such a distribution applies 
for X-ray outbursts in galactic nuclei. Low-mass fueling events, such as partial tidal 
disruptions or the accretion of brown dwarfs, planets, or small gas clouds, would then 
be more common than high-mass fueling events. 

In this paper, we utilize data from the Chandra Deep Fields (see Brandt & Hasinger 2005 for 
a review) to constrain the rate of X-ray outbursts with harder X-ray spectra, higher redshifts, 
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and lower luminosities than those studied to date. Observations of the Chandra Deep Fields 
were made in several well-separated and individually sensitive "epochs," allowing effective 
constraints to be placed upon X-ray outbursts that evolve on timescales of months in the 
rest frame. Some variability work has been performed on the Chandra Deep Field sources. 
For example, Paolillo et al. (2004) studied the X-ray variability of ~ 350 sources detected 
in the Chandra Deep Field-South (CDF-S), but this work has not been optimized to detect 
and to constrain systematically X-ray outbursts of the type relevant here. 

We adopt H = 70 km s _1 Mpc -1 , Qu = 0.3, and = 0.7 throughout this paper. All 
coordinates are J2000. The relevant Galactic column densities are 1.3 x 10 20 cm -2 for the 
Chandra Deep Field-North (CDF-N; Lockman 2004) and 8.8 x 10 19 cm' 2 for the CDF-S 
(Stark et al. 1992). All X-ray fluxes and luminosities quoted throughout this paper have 
been corrected for Galactic absorption using these column densities. 



THE SURVEY 



I n this survey, we focus ed on optically d etected galaxies in the CDF-N (jBrandt et al. 



200ll : Alexander et all 120031 ) and the CDF-S flGiacconi et all 12002 ; Alexander et al 



20031 ). 



An optical galaxy sample was compiled from several catalogs, with the redshifts collected 
from literature. We grouped the Chandra observations into several well-separated "epochs" . 
The five epochs for the CDF-N span 798 days, and the four epochs for the CDF-S span 1 828 
days. At the median redshift of ~ 0.8 for the galaxies in our sample, the rest-frame coverage 
would be 443 days for the CDF-N and 1 016 days for the CDF-S. X-ray count rates or upper 
limits for every optical galaxy were measured for each epoch in three standard Chandra 
bands: 0.5-8.0 keV (full band; FB), 0.5-2.0 keV (soft band; SB), and 2-8 keV (hard band; 
HB). These measurements were then analyzed to search for any X-ray outburst candidates. 
Here we define "X-ray outbursts" to be transient events in galaxy nuclei that cause the count 
rate to vary by a minimum factor of 20 in one of the three standard bands. A minimum 
variability factor of 20 was chosen to discriminate against normal AGN variability, which 



typically has variability factors of 
2004). 



i2— 5 but can be as large as ^10-15 (e.g., Paolillo et al. 



2.1. Optical Galaxy Sample 



Optical galaxies were selected from the catalog of ICapak et al.l (120041 ) for the CDF-N, 
and from the COMBO-17 catalog JWolf et al.l l2004f ) for the CDF-S. To ensure that these 
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galaxies were covered by either the CDF-N or the CDF-S, we used the 0.5-8.0 keV exposure 
maps of the 2 Ms CDF-N and 1 Ms CDF-S as filters. A small portion (~ 30 arcmin 2 ) of 
the CDF-N was not covered by the Capak catalog, so the coverage of this survey is ~418 
arcmin 2 for the CDF-N and ~391 arcmin 2 for the CDF-S. We only chose optical galaxies 
with AB magnitudes R < 25, to optimize our selection of a large galaxy sample with redshift 
measurements. We examined manually those sources with R < 22 and removed some false 
sources, which were generally around bright objects. Galaxies in the COMBO- 17 catalog 
with values for photometry flags > 8 were ignored, as suggested by the creators of that 
catalog. In this way, we selected 13 699 galaxies in the CDF-N and 11077 galaxies in the 
CDF-S. 

We searched the literature for redshift measurements for this op tical galaxy s ample . 
Spectroscopi c redsh ifts were obtained from the spectroscopic surveys of lCowie et al.l (20041) . 
Reddv et all (bood ). and the "Team Keck" Treasury Reds hift Survey Jwirth et all l2004h 



VLT/FORS2 spectroscopic s urvey (|Vanzella et al.l 120051 . 120061 ) in the CDF -S. Photometric 



redshi fts w ere obtained from ICapak et al.l (120071 ) in the CDF-N, and from iMobasher et al. 



in the CDF-N, and from the VI MOS VLT Deep Surv ey flLe Fevre et al.l |200~J" and the 



(2004 


) and 


Wolf et al. 


(2004) 



correlating different optical catalogs. We did not find redshift measurements for 108 faint 
galaxies in the CDF-S. After removing these sources our final optical galaxy sample contained 
13 699 and 10 969 galaxies in the CDF-N and CDF-S, respectively. The spatial distributions 
of these galaxies in the CDF-N and CDF-S are shown in Figure [TJ We found spectroscopic 
redshifts for 2 100 galaxies in the CDF-N and 1 534 galaxies in the CDF-S, with the remaining 
sources having photometric redshifts. The redshift distributions are shown in Figure [2j To 
test the reliability of the photometric redshifts, we examined sources with both photometric 
and spectroscopic redshift measurements, and we studied the offset distribution of these two 
types of redshifts in the CDF-N and CDF-S separately. In both fields, the median offsets 
approach 0. The interquartile ranges are only 0.11 in the CDF-N and 0.20 in the CDF-S. 
The "stellarity index" (class_STAR parameter in SExtract or) is a robust galaxy ide ntifier, 



with for confirmed galaxies, and 1 for confirmed stars (e.g.. lGroenewegen et al.ll2002l ). Less 
than 2.5% of these galaxies had stellarity indices > 0.9, so stars were well removed from this 
sample. 

This sample of 24 668 galaxies contains a representative mix of field g alaxies. Only 



2% of the galaxies are X-ray detected, and thus the AGN fraction is small. iLehmer et al. 



(120071 ) constructed a similar galaxy sample in the Chandra Deep Fields, with z 85Q < 23 
and z = 0-1.4. They studied the rest-frame optical colors as well as the Sersic indices (e.g., 



Haussler et al.l 120071 ) of galaxies in their sample, finding 2 544 late- type galaxies and 727 
early-type galaxies out of 3 271 galaxies. Our sample should have a similar composition 
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in terms of galaxy morphology. Recent close-pair studies have suggested that the merger 

7% 



Bundv et all 12004 . 



rate of field galaxies re mains fairly smal l out t o z ~ 1 (e.g. 
infrared-selected pairs). lElmegreen et al.l (120071 ) examined the GEMS (Galaxy Evolution 
from Morphology and SEDs) and the southern GOODS (Great Observatories Origins Deep 
Survey) fields, which overlap with the CDF-S, and identified ~ 300 interacting galaxies to 
z ~ 1.4 out of over 8 000 optical galaxies. Thus our sample should also be dominated by 
non-interacting galaxies. 



2.2. X-ray Data 



X-ray data analyses wer e performed on the C DF-N, CDF-S, and the Extended Chandra 



Deep Field-South (E-CDF-S; |Le 



11 observations for the CDF-S (Alexander et al. 



imer et al.ll2005). The re are 20 observations for the CDF-N, 



20031 ). and 9 observations for the E-CDF-S 



( iLehmer et al.l 120051 ). We grouped these into "epochs" according to their observation dates, 
five epochs for the CDF-N and four epochs for the CDF-S (see Table [I]). The E-CDF-S 
observations took place ~ 4 years after the last CDF-S observation and were counted as the 
fourth epoch for the CDF-S. The E-CDF-S covered the entire CDF-S and consisted of four 
distinct observational fields with different observation dates. We list these fields separately 
in T able ffl We created images of each epoch from the reduced and cleaned level 2 event 
files ( [Alexander et al.ll2003l ; ILehmer et al.ll2005l ). using the standard ASCA grade set (ASCA 
grades 0, 2, 3, 4, and 6) for the three standard bands. There were thus 15 CDF-N images, 9 
CDF-S images, and 12 E-CDF-S images. The aim point for a given image was taken to be 
the average value over its observations weighted by exposure time. 

To construct X-ray source lists, we ran wavdetect ( Freeman et al.l I2OO2I ) on the 24 
images in the CDF-N and CDF-S using a "a/2 sequence" of wavelet scales (i.e., 1, a/2, 2, 
2 a/2, 4, 4 a/2, and 8 pixels). The false-positive probability threshold in each wavdetect 
run was set to 1 x 10~ 6 . Source lists for the E-CDF-S were taken directly from ILehmer et al. 



(120051 ) because their source-detection method is the same as that used here. To improve 
the accuracy of the X-ray source positions, we matched the optical and X-ray source lists 
with a matching radius of 2'.'5 and centered the distributions of offsets in right ascension and 
declination between the optical and X-ray source positions. This resulted in small (< l'.'O) 
image-dependent astro metric shifts for all X-r ay sources. Compared to the main Chandra 
X-ray source catalogs in I Alexander et al.l (120031 ) . our X-ray source lists give ~ 10 new sources 
for each image, due to the less-conservative false-positive probability threshold we were using 
or the variability of some sources. 



We are interested in the X-ray properties of the optical galaxies selected in §2.1, so 
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Fig. 1. — Spatial distributions of the galaxies in this survey in the CDF-N and CDF-S. Dots 
with different colors represent galaxies in different magnitude bins. Squares are galaxies with 
X-ray counterparts as listed in Table [fj The gap at the top of the CDF-N is due to the lack 
of coverage of the optical catalog. Color figures are available online. 
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Fig. 2. — Redshift distributions of the galaxies in this survey for the CDF-N and CDF-S. 
The median and mean redshifts are 0.78 and 0.91 for the CDF-N, and 0.84 and 1.06 for the 
CDF-S, respectively. The dashed lines show the median values. 
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we searched for X-ray counterparts of these galaxies by cross-correlating the optical and 
X-ray source lists with a matching radius of l'.'O for sources within 6' of the average aim 
point and l'/3 for larger off-axis angles. In 47 cases where there is more than one optical 
galaxy associated with an X-ray source, the closest one is selected. None of the galaxies is 
associated with more than one X-ray source. A summary of these findings is given in Table [TJ 
Around one third of the X-ray sources do not have counterparts in the galaxy sample, either 
because their R magnitude is fainter than 25 or they are stars. We also estimated the 
expected number of false matches by artificially offsetting the X-ray source coordinates in 
right ascension and declination by 5'.'0 (both positive and negative shifts) and re-correlating 
with the optical sources. On average, the number of false matches is only ~ 6% of the total 
matches in our sample for both the CDF-N and CDF-S. 



2.3. Source Characterization and Outburst Searching 

We measured the X-ray counts coincident with every optical galaxy in all 36 images using 
circular-aperture photometry. We have chosen the aperture radii based on the encircled- 
energy function of the Chandra PSF, which was calculated using MARX (Model of AXAF 
Response to X-rayd3) simulations. The 90% encircled-energy radius of the PSF was used. The 
circular aperture was centered at the position of every optical galaxy or its X-ray counterpart 
if it was X-ray detected. If a galaxy was not fully covered by a certain image, then it was 
marked as undetectable for this image. 

The local background was determined in an annulus outside of the source-extraction 
region using background maps with known X-ray sources carefully removed to avoid possible 
contamination from adjacent sources. To create a background map for a given image, we 
first merged our Chan dra X-ray source list from this image to the main Chandra catalogs 



([Alexander et aJj 120031 ) to get a combined source list, using a matching radius of 2'.'5 for 
sources within 6' of the average aim point and 4'/0 for larger off-axis angles. Then we masked 
out these sources using apertures with radii twice that of the 90% PSF encircled-energy 
radius. We filled the masked regions for each source with a local background estimated by 
making a probability distribution of counts using an annulus with inner and outer radii of 
2 and 4 times the 90% PSF encircled-energy radius, respectively. The local background for 
every optical galaxy was then determined using the same annulus on this background map. 

When an optical galaxy was X-ray detected (i.e., had an X-ray counterpart) in a given 
image, the net number of source counts was calculated by subtracting the expected num- 



1 See |http://space.mit.edu/CXC7M ARX 



Table 1. Observation Epochs for the Chandra Deep Fields 







Exposure Time 


Obs. Date a 


Number of Sources b 


Epochs 


Obs. IDs 


(ks) 


tobs 


FB 


SB 


HB 


CDF-N Epoch 1 


580, 967, 966, 957 


221.6 


1999 Dec 12 


131 


105 


73 


CDF-N Epoch 2 


2386, 1671, 2344 


267.4 


2000 Nov 22 


125 


1 r\r\ 

100 


76 


CDF-N Epoch 3 


2232, 2233, 2423, 2234, 2421 


488.2 


2001 Feb 26 


142 


118 


96 


CDF-N Epoch 4 


3293, 3388, 3408, 3389 


385.7 


2001 Nov 17 


138 


115 


82 


CDF-N Epoch 5 


3409, 3294, 3390, 3391 


581.7 


2002 Feb 17 


147 


128 


92 


CDF-S Epoch 1 


1431-0, 1431-1 


114.7 


1999 Nov 15 


77 


68 


48 


CDF-S Epoch 2 


441, 582 


186.5 


2000 Jun 01 


83 


68 


53 


CDF-S Epoch 3 


2406, 2405, 2312, 1672, 2409, 2313, 2239 


638.2 


2000 Dec 18 


136 


124 


81 


CDF-S Epoch 4 F01 c 


5015, 5016 


240.1 


2004 Mar 01 


52 


44 


39 


CDF-S Epoch 4 F02 c 


5017, 5018 


227.4 


2004 May 15 


40 


36 


29 


CDF-S Epoch 4 F03 c 


5019, 5020 


240.7 


2004 Nov 17 


22 


18 


13 


CDF-S Epoch 4 F04 c 


5021, 5022, 6164 


246.0 


2004 Nov 16 


26 


24 


18 



Observations within each epoch only lasted for a few days, so we neglected the length of the epochs and only show 
the average observation dates weighted by exposure time. 

b Number of galaxies in the sample that are X-ray detected in this epoch for the particular band: 0.5-8.0 keV (FB), 
0.5-2.0 keV (SB), or 2-8 keV (HB). 

c Observational helds of the E-CDF-S. Note that a large fraction of each of these helds extends outside the region of 
our outburst survey. 
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ber of background c ounts from the number of counts in the aperture. Poisson errors were 
calculated following iGehrelsl (119861 ) and were propagated through this calculation. When 
an optical galaxy was not detected, an upper limit was calculated. If the number of counts 
in the aperture was < 10 , the u pper limit was calculated using the Bayesian method of 
Kraft. Burrows, fc Nousekl (|199ll ) for 99.87% confidence. For a larger number of counts in 
the aperture, a 3cr upper limit was set by multiplying the square root of the number of 
background counts by three. The number of source counts was then divided by the effective 
exposure time, which was the average value within the aperture on the exposure map, to get 
the X-ray count rate (this procedure corrects the count rate for vignetting and other effects). 

The E-CDF-S, which serves as the fourth epoch for the CDF-S, contains four obser- 
vational fields, which overlap in a few areas over ~50 arcmin 2 . If an optical galaxy was 
detected in more than one field, we chose data from the field with greater source counts. If 
it was detected in only one field, we chose data from this field. If it was not detected in any 
field but was in the overlapping area, we chose data from the field with the longest effective 
exposure time. 

We searched for outbursts following these steps: 



1. For every optical galaxy, we calculated its count rate or upper limit on count rate in 
each epoch and each energy band (i.e., each of the 36 images), unless it was not fully 
covered by a given image. A galaxy had coverage in up to five epochs if it was in the 
CDF-N and up to four epochs if it was in the CDF-S. 

2. For each energy band, we selected galaxies that were detected in at least one epoch. 
For each of these galaxies, we compared its highest count rate with its lowest count 
rate or upper limit on count rate. Thus we got either the variability factor or its lower 
limit for each source. If there was a variation of more than a factor of 20, then we 
considered that there was a candidate for an X-ray outburst in this galaxy. 



3. RESULTS AND DISCUSSIONS 

After systematically analyzing the count-rate variations, we found no outbursts in either 
the CDF-N or in the CDF-S. The distribution of the variability factors in the SB and HB is 
shown in Figure [3j The median relative uncertainty of these factors is ~ 30%. The count- 
rate variations of a few sources exceed a factor of 10. However, these are all off-axis sources 
and the significances of the variations are < lcr. 
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3.1. Upper Limits on the Outburst Event Rate 



The nondetection of any X-ray outbursts in this survey can be used to constrain the 
rate of such outbursts in the Universe. Since the detectability of an outburst depends upon 
its luminosity, the constraints will have a luminosity dependence. Moreover, not all galaxies 
in this survey are capable of producing outbursts. Galaxies without a central SMBH (e.g., 
dwarf irregulars) cannot tida lly disrupt stars, an d a SMBH of mass Mbh ^ 3 x 1O 8 M 
will swallow stars whole (e.g.. lFrank fc Reeslll976l ). SM BH candidates with masses of a few 



1O^M in the centers of galaxies have been found (e.g., iGreene fc Hoi 12004 ; iPeterson et al. 



20051 ). Thus we consider that only galaxies with a SMBH mass greater than 1O 5 M and 
less than 3 x 1O 8 M can produce outbursts. SMBH masses were roughly estimated for all 
galaxies in this survey with the relation between Mbh and total galaxy luminosity in the 
K band, LK,totai- The relation was de rived from the da t a for 27 low-redshift galaxies (10 
late- type and 17 early-type galaxies) in lMarconi fc Huntl (120031 ). The intrinsic dispersion of 
this relation is ~0.5 dex in log Mbh- Although there is evidence showing cosmic evolution 
of the Mbh-0" (S MBH mass and bulge velocity dispersion) relation and the bulge-to-SMBH 
mass ratio (e.g.. lWoo et al.ll2006l ). the coexistent luminosity evolution makes the SMBHs at 
z > 1 coincidently fall on nearly th e same Mbh vers us i?-band magnitude (Mr) relation 
(to 0.3 mag) as low-redshift galaxies (jPeng et al.l 120061 ) . Thus we can expect that the M B h- 
^K,totai relation also holds approximately at z ~ 0.8 for this galaxy sample. L K>totai for 
the galaxies in th is survey was extrapolated from HK' -hand magnitudes (10826 galaxies, 
Capak et al.ll2004h or z'-band magnitudes f287 3 galaxies. ICapak et al.ll2004r) for the CDF-N, 
and from z-b and magnitudes (10167 galaxies, ICaldwell et al.l 120051 ) or _R-band magnitudes 
(802 gala xies, Wolf et all 2004) for the CDF -S, using the spectral energy distribution for Sbc 
galaxies (jColeman. Wu. &: Weedmanlll980l ). The resulting conversions from HK', z', z and 
R magnitudes (AB system) to the standard K magnitude were given by HK 1 — K = 1.65, 
z' — K = 2.07, z — K = 2.01, and R — K = 2.70. Based on the dispersion of the MBH--^K,totai 
relation and the uncertainties in the color conversions, the derived Mbh is estimated to be 
good within an order of magnitude (as we show below, our main results are not sensitive to 
Mbh)- A comparison between the SMBH mass functions derived for the low-redshift galaxies 
in this sample and the local SMBH mass function for all galaxy types (jMarconi et al.l 12004 ) 
is shown in Figure HI There is basic agreement between the shapes of these mass functions 
indicating that our mass measurements are reasonable. 

There are 11339 galaxies in the CDF-N and 8 268 galaxies in the CDF-S with Mbh hi 
the range of 1O 5 M < Mbh < 3 x 1O 8 M . There will also be another constraint set by the 
physics of outburst production, i.e., a SMBH with a given mass cannot produce outbursts 
with arbitrarily high luminosities owing to, e.g., the Eddington limit. However, we ignore 
this constraint for now, as we would first like to present model-independent limits. 
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Variability Factor Variability Factor 

Fig. 3. — Histograms showing the distribution of variability factors in the SB (upper panels) 
or HB (lower panels) for galaxies in the CDF-N and CDF-S. The variability factor of a given 
galaxy is defined as the ratio between its highest count rate and its lowest count rate or 
upper limit on count rate. The shaded areas are for galaxies with lower limits on variability 
factors only. 
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Fig. 4. — SMBH mass functions f or galaxies with redshift less than (a) 0.5 or (b) 1.0, and the 
local SMBH mass function from iMarconi et al.l (120041 . solid curve) scaled to compare with 
the mass functions for our sample. The dotted curves indicate mass functions in the CDF-N 
while the dashed curves are for the CDF-S. The deviation at low SMBH mass in the right 
panel is simply due to a selection effect; it is difficult to detect the galaxies hosting SMBHs 
with small masses at relatively high redshift s. 
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To determine the allowed rate of outbursts of X-ray luminosity greater than or equal to 
a certain value, we used the following "recipe": 

1. For a given energy band (FB, SB, or HB), we picked an X-ray luminosity, Lx,burst> as the 
lower limit, and we derived the rate of outbursts of X-ray luminosities > Lx.burst- Count 
rates or upper limits on count rates for all the galaxies in the survey were converted 
to luminosities based on their redshifts, using the Portable, Interactive, Multi-Mission 
Simulator (PIMMS) for a power-law photon index of T. Galaxies generally have soft 
X-ray spectra during outburst and harder spectra in the their basal states. We thus 
considered a few choices of T ranging from 2 to 5. The luminosity distributions for 
T = 2 and T = 4 are shown in Figure Most of the luminosities are upper limits. 

2. We assumed that the typical duration of an outburst was T burst . Theoretically predicted 
light curves of outbursts show a characteristic fast rise and slow decay. After the bulk of 
the material is accreted, which could be on a timescale of the order of months, the debris 
starts to fo rm a radiation-pr essure supported torus and the luminosity declines slowly 



as oc t 5//3 ( Reeslll988l . Il990l ). This long-term decay could last for ye ars and may have 



been observed (with large light-curve gaps) in a few outbursts (e.g.. iKomossa &: Bade 



19991 ). The detailed luminosity evolution of outbursts, especially for the first few 
months, remains unknown, and likely varies greatly from event-to-event. For simplicity, 
we assumed that the luminosity is constant during the outburst, and after that the 
X-ray luminosity would drop by a minimum factor of 20. We then determined the total 
rest-frame time over which we are sensitive to outbursts of luminosity > -Lx,burst : 

n 

Ttotal ,burst> ^burst) — ^ Tj scns (Lx, burst ; ^burst) ; (1) 

i=l 

where n = 19607 is the number of galaxies in the survey with 10 5 M Q < Mbh < 3 x 10 8 M q , 
and Tj )Sens is the rest-frame time over which we are sensitive to an outburst of luminosity 
> ^x.burst for galaxy i. T i>sens can be expressed as: 



Ti sens — / Wi(t)dt. (2) 

*/Time,rf 

This integration is over all the rest-frame observation period, from time Tb urs t before 
the first epoch to the last epoch, and Wi(t) is our sensitivity "window function" to 
outbursts starting from time t and with luminosity > Lx,burst and duration T burst . If 
we could detect such an outburst based on our searching criteria (step 2 in §2.3), then 
Wi(t) = 1; otherwise Wi(t) = 0. 
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3. The 90% confidence upper limit on events is 2.303 from iGehreld ( 119861 ). Thus the 



90% confidence upper limit on the event rate of outbursts is given by 

• 2.303 outbursts 

-/total galaxy yr 

4. Repeat the above steps for various values of Lx,burst and T burst as well as different 
assumptions about outburst spectral shape. 

The constraints set by Equation [3] depend significantly on the outburst duration T burst 
and the lower luminosity limit Lx,burst, and they weakly depend on the spectral shape (we 
consider X-ray photon indices T =2, 3, 4, or 5) and energy band (the SB or HB), as shown 
in Figures [6], [7| and [HJ Longer outburst durations lead to tighter upper limits on the event 
rate. Photon indices affect the derived X-ray luminosities through the conversion from 
count rate to flux (from PIMMS) and the K correction term (1 + z) <yT ~ 2 \ As V increases, the 
conversion factor increases in the SB and decreases in the HB, while the K correction term 
always increases and preferentially dominates over the conversion factor for high-luminosity 
sources since they generally have large redshifts. Thus the dependence of the upper limits 
on the spectral shape behaves differently in the SB and HB. Note that luminosities Lx,burst 
in the SB and HB are different, so the upper limits on the event rate in these energy bands 
cannot be compared to each other directly. Generally, assuming an outburst duration of 
6 months, which was also adopted in Donley et al. (2002), the upper limit on the event rate 
is ~ 10~ 4 galaxy -1 yr -1 , for an outburst with X-ray luminosity > 10 43 ergs s _1 . 



3.2. Comparison with Previous Results 



We compared these new constraints with the theoretical study by IWang &: Merritt 



( 120041 ) using a singular isothermal sphere, which analytically predicted a stellar tidal-disruption 
rate of 



7 /2 / A/ r \ -1 / m \ -1/3 / p x 1/4 



(4) 

with Mbh being the SMBH mass, a the velocity dispersion of the host galaxy, and m* 
and R+ the mass and radius of the tidally disrupted stars. We added an "amplification" 
factor A(z) here to represent any redshift evolution of the rate. As the evolution of X-ray 
outbursts is unknown (s ee §1), we set A(z) = 1 for now. Utilizing the Mbh-c relation from 



Ferrarese fc Fordl (120051 ): 



M BH = 1.66x 1 0«M e ( 5sn ^- =T ] (5) 



4.86 
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Fig. 5. — SB and HB X-ray luminosity distributions of galaxies in the CDF-N and CDF-S. 
Luminosities for all the epochs were collected, so each galaxy is plotted multiple times, 
corresponding to the number of epochs in which it was observed. Most of the luminosities 
are 3a upper limits. Luminosities of X-ray detected galaxies are represented by the shaded 
area. The number of galaxies was plotted using a logarithmic scale to show the small fraction 
of X-ray detected galaxies. A photon index of T = 2 or T = 4 was adopted in these plots. 
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Fig. 6. — Dependence of T total on outburst duration at different outburst luminosities. T total 
is the total rest-frame time over which we are sensitive to outbursts of (a) SB luminosity 
-^SB.burst or (b) HB luminosity Lhb, burst; see Equation CD A photon index of T = 2 was 
adopted when making these plots. 
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Fig. 7. — The same as Figure El but for a photon index of T = 4. 
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Equation @] becomes 

/ », \ -0.28 / \ -1/3 / D \ 1/4 

As we are here comparing with a physical model, we modified our previous constraints 
on event rate by also considering the capability of a SMBH with a given mass to pro- 
duce outbursts of high luminosities. We simply assumed that there was a narrow range 
for the outburst luminosity, which was around a fraction / B dd of the Eddington luminos- 
ity /Edd = -^boi/^Edd, where Lboi is the outburst's bolometric luminosity and LEdd ~ 
1.3 x 1O 38 (Mbh/M ) ergs s -1 is the Eddington luminosity. We also assumed a bolometric 
correction f hc , which was defined as /be = A>oi/-kx,burst- The constraint on M B h is then 

Mbh > J hcL ^T M 0- (7) 

1.3 x 10 38 /Edd 

Applying this SMBH-mass requirement in addition to the 1O 5 M < Mbh < 3 x 10 8 M q 
requirement in step 2 of §3.1, we got weaker constraints on the event rate for high outburst 
luminosities, since the number of qualified galaxies n is smaller. The derived upper limits 
as a function of Lx,burst for the SB and the HB are shown in Figures M and [101 respectively. 
Here we assume T burst = 6 months, /be = 10, and f^dd = 1-0 or 0.1. The event rate shows a 
dependence on the Eddington ratio, /Edd, since smaller / E dd will limit the number of galaxies 
that are capable of making bright outbursts. The new requirement on the SMBH mass does 
not affect the outburst rate at small Lx,burst, because all galaxies are capable of producing 
outbursts of such low X-ray luminosity. 

To employ the formula for the theoretical tidal-disruption rate of a single galaxy (Equa- 
tion ED in our survey, we used Mbh calculated in §3.1. As the dependence of Nwu on Mbh 
is weak, there should not be a large error if our SMBH-mass estimation is not highly ac- 
curate; an error of a factor of 5 in Mbh only changes A^vm by a factor of 1.5. As we are 
only considering tidal disruptions of stars her e, the disruption rate is d ominated by subsolar 



stars according to the stellar mass function (jMilosavlievic et al. 



2006h. The radius of these 



stars follows the relation R*/R & « (m*/M ) - 8 (e.g.. iKippenhahn fc Weigertlll990l ). Thus 
the dependence on stellar mass and radius is weak, and we simply assumed solar mass and 
radius. For each Lx,burst, we took the average value of N WM f° r & U galaxies for which we are 
sensitive to an outburst of luminosity > L Xi burst, which required T ijSens > 0. The predicted 
event rate for this survey as a function of L x ,burst is then given by 

• EILi^wm outbursts 

2^=i k i galaxy yr 
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Fig. 8. — 90% confidence upper limits on the event rates of outbursts derived from this 
survey as functions of the lower limit on (a) SB or (b) HB X-ray luminosity; see Equation 
El Solid, dotted, dashed, and dash-dotted lines represent photon indices T of 2, 3, 4, and 5, 
respectively. We show the constraints on two different outburst durations, 1 month and 6 
months. We made no assumptions about the physical process causing the outburst except 
that the SMBH mass is in the range 1O 5 M < M BH < 3 x 1O 8 M . 




Fig. 9. — 90% confidence upper limits on the event rates of outbursts derived from this 
survey and expected event rates from the theoretical prediction as functions of the lower 
limit on the SB X-ray luminosity of outbursts. Solid, dotted, dashed, and dash-dotted lines 
represent photon indices T of 2, 3, 4, and 5, respectively. SMBHs that were not capable 
of producing outbursts of luminosity LsB,burst were removed from the sample. We assumed 
^burst = 6 months, /be = 10 and /ecM = 1-0 (a) or 0.1 (b). The predicted rates were derived 
under the assumption of isothermal stellar density distributions in galactic nuclei. 
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where hi — 1 if T ijSens > 0, and fcj = if T,- 



i, sens 



0. 



iV t h as a function of Lx.burst for different spectral shapes and Eddington ratios is also 
plotted in Figures M arid [TUJ Comparing these predictions with our survey constraints, we 
see that the upper limits derived in this survey are consistent with the predicted rates, 
except that around an X-ray luminosity of 10 43 ergs s _1 , this deep- field survey sets a tighter 
constraint on the rate of Eddington-limited events than the theory (Figure [9^, Figure ITUa). 
If we reduce /Edd to 0.1, the discrepancy decreases and almost disappears (Figure 03, Figure 
[TUb). There are other sources that may help to resolve this discrepancy: (1) As mentioned in 
§1, a tidal-disruption event may not exhibit the expected outburst characteristics, due to a 
short duration of the event or X-ray obscuration. It is then likely to be missed by current and 
previous surveys. (2) Equation @] is derived under the assumption that the galactic nucleus is 
a singular isothermal sphere. For other k inds of density distribut ions, the predicted rate will 
become smaller according to Figure 5b of I Wang fc Merrittl (120041 ). which gives the computed 
tidal- disrupt ion rates and SMBH masses for all the galaxies in their sample. A straight-line 
fit to the data points in this plot gives 



WM-Rev 



BH, 



2.3 x 10~ 4 yr 



-i 



Mi 



(l0 6 M 6 



-0.52 



(9) 



with a smaller rate and stronger dependence on SMBH mass. A comparison with the pre- 
dicted tidal-disruption rate derived from this modified relation is shown in Figure HU The 
revised rates iVth-Rev are ~ 0.5 dex smaller and now are slightly below our observed upper 
limits for the full luminosity range. (3) The Mbh-c relation in Equation [5] only hol ds for local 



SMB Hs. The relation could evolve in the sense of velocity dispersion decreasing (IWoo et al. 



20061 ) with redshift for a fixed BH mass, making the predicted rate decrease at high red- 
shift. All these effects may cause the upper limits on event rate given by this survey to be 
above the predic t ed ste llar tidal-disruption rates. On the other hand, the galaxy sample in 
Wang fc Merrittl (120041 ) is a set of 61 elliptical galaxies, while our sample is dominated by 
late-type galaxies. By considering the contribution from the b ulges of spirals, the predicted 
rates could be significantly increased ([Wang fc Merrittl I2004J ). as the masses of the bulges 
are smaller than the total masses of the galaxies that we were using. This could increase the 
discrepancy between the observations and theoretical predictions. 

We employed the Eddington limit above based on the assumption that the X-rays come 
from a transient accretion disk and its corona. Alternatively, X-rays could instead come 



from tidal-stream collisions (e.g.. lKochaneklll994l ). in which case the luminosity could per- 



haps exceed the Eddington limit. If we remove the constraint on the SMBH mass set by 
Equation [7J the upper limits on the event rates get tighter for high X-ray luminosities and 
are the same as those in Figure [BJ and the predicted rates also change due to their depen- 
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Fig. 10. — The same as Figure EE but for the HB X-ray luminosity of outbursts. 




41 48 43 44 45 41 42 43 44 45 

lOg L SB,b„r«t. ( er S S S ~') l°g L HB,bur»t ( er g S ^) 



Fig. 11. — 90% confidence upper limits on the event rates of outbursts derived from this 
survey and expected event rates from the theoretical prediction for the (a) SB and (b) HB, 
under the assumptions of T burst = 6 months, /be = 10 and /Edd = 1-0. Solid, dotted, 
dashed, and dash-dotted lines represent photon indices T of 2, 3, 4, and 5, respectively. The 
observational upper limits are the same as those in Figures [9] and [101 while the predicted 
rates were derived from a modified relation which does not require isothermal stellar density 
distributions in galactic nuclei (Equation [9]). 
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dence on Mbh- The results are plotted in Figure [T2| which shows the predicted rates from 
both the original relation (Equation [6]) and its revision (Equation [9]). Our survey constraints 
for Lx.burst ^ 10 43 ergs s -1 are ~ 0.5 dex tighter than the analytical predictions assuming 
isothermal stellar density distributions in galactic nuclei, an d are consistent w i th the revised 



tidal- disrupt ion rates derived from the full galaxy sample in IWang fc Merrittl (120041 ) . 



We have assumed A(z) = 1 in the above analyses, which represents no redshift evolution 
of the outburst rate. However, the rate of X-ray o utbursts could increase wit h redshift owing 



to changes in galactic nuclei and SMBH masses. iMilosavljevic et al.l (120061 ) have proposed 
that stellar tidal disruptions are largely responsible for the AGN X-ray luminosity function at 
luminosities below 10 43 — 10 44 erg s _1 . If this is indeed the case, then the order-of- magnitude 
increase in the comoving number density of such AGNs out to z ~ 1 (e.g., Brandt & Hasinger 
2005 and references therein) implies that the outburst rate must correspondingly increase. 
Scaling the Wang & Merritt (2004) rates upward by A(z = 0.8) ~ 10 to account for the 
median redshift of our sample would lead to significant disagreement with our observational 
constraints, even if the revised rate (Equation [9]) is adopted, suggesting that the X-ray 
luminosity function at luminosities below 10 43 — 10 44 erg s _1 is not primarily due to stellar 
tidal disruptions. 

On the observational side, Donley et al. (2002) performed a systematic survey for X-ray 
outbursts using the ROSAT database (although this survey had some substantial systematic 
uncertainties owing to complex selection effects). They detected five outbursts and placed 
the first constraints on the rate of such outbursts. The rate of large-amplitude X-ray out- 
bursts from inactive galaxies in the local Universe is ~ 10~ 5 galaxy^ 1 yr _1 , estimated from 
the survey volume and the galaxy space density. Compared to their results, our survey con- 
straints are based on uniform observational data for 24 668 optical galaxies which are less 
biased, and there are no uncertainties introduced in the estimation of survey volume or the 
galaxy space density. Moreover, we are able to probe X-ray outbursts with higher redshifts 
and in a harder X-ray band for the first time, and we derived luminosity-dependent rate 
constraints which offer insight into the low-luminosity regime. HB X-ray surveys are able 
to detect obscured outbursts which could have been missed in previous surveys; intrinsic 
column densities as low as ~ 5 x 10 21 cm -2 , perhaps associated with gas from the tidally dis- 
rupted star, would greatly reduce the detectability of soft-spectrum outbursts in the ROSAT 
band. Here we adopted Tb urst = 6 months, the same as in Donley et al. (2002). The SB 
and HB constraints are shown in Figure [HJ From Lx.bm-st ~ 10 43 down to 10 41 ergs s _1 , the 
upper limit on the event rate increases from ~ 10~ 4 galaxy -1 yr _1 to ~ 10~ 2 galaxy -1 yr _1 , 
mainly due to the limited sensitivity in the low-luminosity regime. Further observations 
are required to assess whether low-amplitude events, which could come from partial tidal 
disruptions or the accretion of brown dwarfs, planets, or small gas clouds, are truly more 
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Fig. 12. — 90% confidence upper limits on the event rates of outbursts derived from this 
survey and expected event rates from the theoretical prediction for the (a) SB and (6) HB, 
under the assumption of T burst = 6 months. Solid, dotted, dashed, and dash-dotted lines 
represent photon indices T of 2, 3, 4, and 5, respectively. The X-ray luminosities were not 
constrained by the Eddington limit in these plots. The observational upper limits are the 
same as those in Figure The predicted rates were derived from Equation [6] which is the 
analytic relation for galactic nuclei with isothermal stellar density distributions, and from 
Equation which does not require isothermal stellar density distributions. 
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common than high-amplitude events, as for other transient phenomena in nature. When 
-^x.burst ^ 10 43 ergs s -1 , the upper limit on the event rate is ~ 10~ 4 galaxy -1 yr _1 , either in 
the SB or the HB. Compared to the results in Donley et al. (2002), these constraints allow 
an amplification of the outburst rate up to a factor of 10 at most when considering obscured 
X-ray outbursts and redshift evolution from z ~ to z ~ 0.8. 



4. CONCLUSIONS AND FUTURE WORK 

In summary, we constructed a sample of 24 668 optical galaxies in the CDF-N and 
CDF-S with redshifts obtained from the literature; the median redshift is ~ 0.8. We analyzed 
exceptionally sensitive Chandra observations of these galaxies, which span 798 days for the 
CDF-N and 1 828 days for the CDF-S. We searched for X-ray outbursts with the criterion 
that the count rate varies by a minimum factor of 20 in one of three standard bands. No 
outbursts were found, and thus we set upper limits on the rate of such events in the Universe, 
which depend on the X-ray luminosity of outbursts. If we only consider those galaxies 
hosting SMBHs, and those with SMBHs not massive enough to swallow a whole star without 
disruption, we derive an upper limit on the rate of an outburst with Lx,burst > 10 43 ergs s _1 
and Tt, urst = 6 months to be ~ 10 -4 galaxy" 1 yr _1 (Figure EJ), without any other assumptions 
about the physical model producing X-ray outbursts. Compared to the survey by Donley 
et al. (2002), our survey probes both higher redshifts and harder X-ray energies. The 
outburst rate may increase by a maximum factor of 10 when taking into account both 
obscured X-ray outbursts and redshift evolution from z ~ to z ~ 0.8. We are able to set 
constraints on low- luminosity events down to 10 41 ergs s _1 , though at this luminosity the 
event rate is limited by sensitivity. 



We also compared our constraints to the predicted tidal-disruption rates of lWang fc Merritt 



(120041 ). exploring several possibilities in physical parameter space. If the outburst luminosity 
is limited by the Eddington luminosity with / Edd < 0.1, or the predicted rate is computed 
from a modified relation which does not require isothermal stellar density distributions in 
galactic nuclei, our results are roughly consistent with theoretical predictions. Otherwise, our 
constraints are tighter than the predictions. Moreover, if stellar tidal disruptions are largely 
responsible for t he AGN X-ray lum i nosity function at luminosities below 10 43 — 10 44 erg s _1 



as proposed by iMilosavljevic et al.l (120061 ) , the predicted rate should be scaled upward by 
A(z = 0.8) ~ 10 to account for the median redshift of our sample. This scaling leads to signif- 
icant discrepancies with our observational constraints, suggesting that the X-ray luminosity 
function at low luminosities is not likely to be dominated by stellar tidal disruptions. 



The constraints set by this study could be significantly improved by further deep-field 
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surveys and new missions. For example, the 1 Ms of additional CDF-S exposure starting 
in 2007 September will increase the monitoring baseline significantly and provide ~ 13 new 
observations which could be grouped into several additional epochs for outburst searching 
and other variability studies. Future missions such as Lobster and eROSITA in the soft 
X-ray band and the BHFP in the hard X-ray band will also b e capable of det ecting and 



studying outbursts, benefiting from their large sky coverage. iGrindlayl (12004 ) predicted 
that the Energetic X-ray Imaging Survey Telescope (EXIST) implementation of the BHFP 
w ould detect X-ray outbu rsts out to ~ 100 Mpc at a rate of ~ 30 yr -1 , adopting the rates 
of IWang fc Merrittl (120041 ) . Such missions will put better constraints on the outburst rate, 
the fraction of obscured outbursts, and redshift evolution, enhancing our knowledge about 
the nature of X-ray outbursts. 
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